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AstroFit is an interface adding astrophysical components to programs for fit- 
ting physics beyond the Standard Model (BSM) to experimental data from 
collider searches. The project aims at combining a wide range of experimental 
results from indirect, direct and collider serarches for Dark Matter (DM) and 
confronting it with theoretical expectations in various DM models. Here, we 
introduce AstroFit and discuss first results. 
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1. Introduction 

Various different experiments explore the properties of Dark Matter (DM) 
while a plethora of theories offer explanations to its nature. It is the goal of 
the AstroFit project to constrain DM models by combining experimental 
data from both astrophysics and collider physics, and thus find the best fit 
regions for the parameter space of the underlying theory. AstroFit itself is a 
Fortran program, serving as an interface between programs used in particle 
physics for fitting physics beyond the Standard Model, such as Fittino, 1-3 
and programs like DarkSUSY 4 ' 5 that can be used to calculate theoretical 
predictions for direct and indirect detection experiments. An overview of 
experimental input usable with AstroFit will be given as well as an example 
of how to use AstroFit in combination with the Fittino program. Here, the 
first results from a Constrained Minimal Supersymmetric Standard Model 
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(CMSSM) fit including information from latest collider (i.e. Large Hadron 
Collider (LHC)), direct detection and indirect detection instruments are 
presented. 



2. Complementarity of Experiments 

Indirect searches for Dark Matter concentrate on finding products from DM 
annihilation, such as photons of different energy ranges (e.g. radio, X-ray, 
7-ray), antiprotons, positrons and neutrinos with specially designed exper- 
iments, ranging from ground-based Cherenkov telescopes like H.E.S.S., 6 
MAGIC 7 and VERITAS 8 over satellite experiments like Fcrmi-LAT 9 and 
PAMELA 10 to balloon experiments as ATIC, 11 among many others. For 
this analysis, we used photon flux upper limits from dwarf galaxies by 
H.E.S.S. 12 and Fermi-LAT as observables, which have been calculated by 
AstroFit as 



d$(Aft,£ 7 ) _ I {av) d7V 7 
AE^ ~ 8^ ml d£ 7 



J(Af2)Af7, (I) 



where d<&(AQ, E 1 )/dE 1 is the differential photon flux, E 1 the photon en- 
ergy, (av) the velocity weighted annihilation cross-section and m x the mass 
of the DM particle, d-/V 7 /d-E 7 the differential number of gammas produced 
per annihilation per energy and J(A£!) = (1/AS1) f An dfl f ; dl p 2 3M (/) 
the integral over the DM density squared along the line of sight. Instead 
of using the photon flux upper limit as direct observable, the thermally 
averaged cross-section times the relative velocity can be used if desired, as 
some experimental limits are preferably presented in terms of (crw) max . 

Direct detection instruments such as DAM A/LIBRA, 13 CoGeNT 14 and 
Xenon, 15 located in underground laboratories, measure signals from DM 
interactions with target elements like Xenon, Germanium or compounds 
such as Nal. Signals are detected via scintillation, phonons or ionization, 
depending on the experiment. Assuming a scattering of a weakly interact- 
ing massive particle (WIMP) as DM candidate with the target material, 
the spin-independent scattering cross-section per nucleon, conventionally 
adopted for comparison between experimental results, is calculated as fol- 
lows in AstroFit: 

^nucleon ' ^ ' 

with a p and a n being the spin-independent cross-section for one proton 
or neutron, respectively, and Z and A being the atomic and mass number 
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of the target element. Latest results from direct detection have recently 
shown a conflict between measurements from different experiments. While 
the DAMA, CoGeNT and CRESST 16 collaborations each published de- 
tections of signals, the Xenon collaboration has shown upper limits lower 
than the signal regions of the previous experiments for the WIMP-nucleon 
cross-section. Possible reasons for this diversity are still being discussed. 

As particle physics input from collider experiments, information on B- 
and Z-physics, like masses, edges in mass spectra, widths, asymmetries, etc. 
have been used in this analysis as well as the anomalous magnetic moment 
of the muon and latest event rates from the ATLAS experiment at LHC 
(see 1-3 for details). Also, the relic density of Cold Dark Matter (CDM) 
provided by WMAP, 17 ft DM /i 2 = 0.1123 ± 0.0035 which can be directly 
compared to the theoretical prediction for thermally produced DM from 
DarkSUSY, is used as an observable in AstroFit. 

3. Structure of AstroFit 

As depicted in figure 1, AstroFit provides an extensive database of rele- 
vant experimental data to add to a fit process, at the moment implemented 
in Fittino. While the steering of AstroFit is done by a user-friendly text 
input file, information on the particle spectrum is given via a SUSY Les 
Houches Accord (SLHA) file directly from Fittino. From this spectrum 
file, AstroFit subroutines are designed to calculate the model predictions 
by using various DarkSUSY functions. The theoretically estimated observ- 
ables are then compared to actual measurements from astrophysical ex- 
periments. From the comparison, a Ax 2 -contribution is calculated for each 
individual observable and is handed to Fittino in each step of the mini- 
mization process to be used together with the information from particle 
physics. In a stand-alone subroutine of AstroFit, the A% 2 is calculated by 
applying A% 2 = Yl ((Ooxp - O t hco)/c OX p) f° r data points or continuously 
in parabolic increase by extrapolation using the given confidence level for 
limits or containment regions for claimed signals for any observable. A min- 
imization strategy in Fittino then determines the lowest global \ 2 1 an d the 
results comprise the best fit regions for all model parameters. 

4. Analysis and Results for CMSSM Model 

Exemplarily, this analysis focussed on the CMSSM, using a Markov Chain 
Monte Carlo algorithm to fit the following parameters (defined at the Grand 
Unification Theory scale): M - the universal scalar mass, Mi/ 2 - the uni- 
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Fig. 1. AstroFit flowchart: Minimization process taking place in Fittino, calculations 
of theoretical observables using subroutines from DarkSUSY and comparison to astro- 
physical observables done in AstroFit. 



versal gaugino mass, Aq - the common trilinear coupling and tan j3 - the 
ratio of the vacuum expectation values of the two Higgs fields. We assume 
the neutralino to be the DM particle and set \i = +1. 

For our analysis, the data from CoGeNT has been used in one fit and 
the data from XenonlOO and Xenongoal in two others (compare 18) . Ad- 
ditionally, fits with data for XenonlT and CRESST 19 are in preparation. 

While claimed DM signals from the CoGeNT experiment could not find 
agreement with a CMSSM fit, upper limits from the Xenon experiment 
give constraints in this scenario. In figure 2, fit results are shown using 
only particle physics information (including latest 2fb _1 results from LHC) 
and the relic density of CDM. Adding further constraints from indirect and 
direct detection, i.e. photon flux upper limits and results from XenonlOO 
data, leads to the results shown in figure 3. In the fit, the bulk region 
is already excluded by all other input conditions, while adding XenonlOO 
limits provides an additional constraint on both the coannihilation and the 
funnel region. 

Photon flux upper limits do not constrain the CMSSM parameter space 
so far. However, using the latest joint likelihood results for dwarf spheroidal 
galaxies from Fermi-LAT (see 20) could show first effects. 



5. Conclusion and Outlook 

Using all available results from DM searches can help confirm, constrain 
or exclude regions in parameter space of DM models remarkably, making 
it possible to edge closer to understanding physics beyond the Standard 
Model in general and the nature of DM in particular. In this study, observ- 
ables from particle physics have been combined with the relic density of 
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Fig. 2. Mo - 2D 2a contour region, 

for a fittino fit using observablcs from col- 
lider production and the relic density of 
CDM. 



Fig. 3. Same as Fig. 2, including also 
XenonlOO lid limits. The additionally 
constrained models lie mostly in the coan- 
nihilation and funnel region. 



CDM, photon flux upper limits from dwarf spheroidal galaxies as well as 
direct detection signals and upper limits, showing considerable impact in 
constraining the CMSSM parameters. While claimed DM signals from the 
CoGeNT experiment could not be fit within a CMSSM model, XenonlOO 
upper limits contributions have been shown. With upcoming fits using As- 
troFit and Fittino, the parameter space of DM models can be constrained 
even further, using latest results from the XenonlOO (as well as predictions 
from Xenongoal and XenonlT) and CRESST direct detection instruments, 
incoming results from the LHC and indirect detection information. As nu- 
merous results in the field of experimental DM physics are expected in the 
near and mid-term future, it is extremely important to strengthen integra- 
tive approaches, which AstroFit facilitates. Within AstroFit, it is planned 
to provide subroutines for all major observables in DM research, such as 
photon fluxes in different energy regimes and from different sources, antipro- 
ton, positron and neutrino fluxes in addition to the relic density of Cold 
Dark Matter and data on the scattering cross-section from direct detection, 
which can be used together with measurements from collider searches. The 
inclusion of available information in particle physics, astroparticle physics 
and cosmology in a combined DM search is an important tool to help inter- 
prete state-of-the-art physics in these disciplines. As such, it is planned to 
release a public version of AstroFit for a larger community of researchers 
to use. 
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